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We show here that the model Hamiltonian used in our paper for ion vibrating in a g-analog 
harmonic oscillator trap and interacting with a classical single-mode light field is indeed obtained 
by replacing the usual bosonic creation and annihilation operators of the harmonic trap model by 
their g-deformed counterparts. The approximations made in our paper amount to using for the 
ion-laser interaction in a (/-analog harmonic oscillator trap, the operator F q = e - ^ -5- ) 'e ltA e teA 
, which is analogous to the corresponding operator for ion in a harmonic oscillator trap that is 

e I 2 . 

F = e~ ( a V ea e lea . Here a and a' are the usual creation and annihilation operators, whereas A 
and A' are the q-deformed bosonic creation and annihilation operators. There is no problem with 
diagonalising this operator using the basis states \g, m) and |e, m) , where m stands for the motional 
number state. In our article we do not claim to have diagonalised the operator F q = e le ^ At+A - ) for 
which the basis states \g,m) and |e,m) are not analytic vectors. 



fiSlIn our article [jlj we have used only q-deformed bosonic creation and annihilation operators, which in the limit 
q — > 1 coincide with the harmonic oscillator creation and annihilation operators. In this reply, in order to clarify 
things we shall make a distinction between the two types of operators and discuss the ion-laser system Hamiltonian 
separately for these two cases. 



A. The Hamiltonian for ion in a harmonic Oscillator trap 

The system that consists of the ion moving in a harmonic oscillator trap potential and interacting with a classical 
single-mode light field of frequency is described by the Hamiltonian, 

H = ~fiij(ala + aa r ) + h\Acr z + ^hn(Fa+ + F^a~) (1) 

where A = uj a — uii , is the detuning parameter and fl is the Rabi frequency of the system. The operator F stands 
for Q| exp(ikx) = exp[ie(a^ + a)]. The parameter e = is a function of the ratio of the recoil energy of the ion 

E r = 2 m anc ^ ^ ne characteristic trap quantum energy E t — huj. Here k is the wave vector of the light field and x 
position quadrature of the center of mass. The second term in the Hamiltonian refers to the energy associated with 



1 



internal degrees of freedom of the ion, whereas the third term is the interaction of the ion with the light field. For 
harmonic oscillator trap one has 



,aU = l; 7V a t _ a l N = a t • Na _ a N 



By using a special form of Baker-Hausdorff Theorem, which reads [[| , 

e X+Y = e -l[X,Y] e X e Y 



(2) 



(3) 



and is valid whenever both X and Y commute with the commutator Z — [X, Y] , the operator F may be rewritten as 
a product of operators i.e 



,(MIV..]) ^ ta (=&) 



F = cxp[ie(a^ + a)] = e v 
Here we examine in some detail as to what kind of processes are implied by the operator 



F 



exp 



(ie) n at™ {ief a k 



fc! 



(4) 



(•5) 



The terms for n > k correspond to the increase in energy by (n — k) quanta while the one's with n < k can destroy 
(k — n) quanta thus altering the amount of energy linked with the center of mass motion. For [n = fc), we have 
diagonal contributions. The contribution from a particular term with operators a* n a k \s determined by the coefficient 

V 2 / V nlk - 



B. The Hamiltonian for ion in a g-analog harmonic oscillator trap 



Next we consider a single two level ion having ionic transition frequency u a in a quantized g-analog quantum 
harmonic oscillator trap(g-deformed harmonic oscillator trap) interacting with a single mode travelling light field. 
The creation and annihilation operators for the trap quanta satisfy the following q-commutation relations, 

AA* - qA^A = q- N ; NA* - A*N = A* ; N A - AN = -A (6) 

Here N is the number operator. The operators A and A* act in a Hilbert space with basis vectors |n), n = 0, 1, 2, 
given by, 

(At)™ 

I") = tttttt 10) (7) 

such that N\n) — n\n) . The vacuum state is A |0) = 0. We define here [x] q as 

(8) 



q - <r 



and the g-analog factorial [n] q \ is recursively defined by 

[0]q! = [l] g ! = 1 and [n] q \ = [n]q[n — l] g l. It is easily verified that 

At | n > = [n + 1]J \n + 1) ; A \n) = [n]\ \n - 1) 

and N is not equal to A + A. 

The g-deformed operators A and A^ are related to a and at through 



(9) 
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A = af{N) ■ A^ = f(N)a^ (10) 
where f(N) = O^k^J 2 and N = a^a. 

The Hamiltonian for an ion interacting with light in a q-analog harmonic oscillator trap may now be written as 

H g = -fiu(AU + AA^) + ^hA(T z + hiVL{F q a+ + F^a~) (11) 



where by analogy with Eq. (0) we choose 



F 9 = eV ' / e ieA V eA . (12) 



which in the limit <? — > 1 reduces to F and can be expanded as 



(-r-) ^ (ie)Mt" ~ (»e) fc A fc 



^ ' 'L,^^^- (13) 



n=Q fc=0 



Various terms in the expansion of this operator represent processes which might result in transitions of the center of 
mass from a given motional state, in the q-analog trap, to another, while loosing or gaining energy. The coefficient 
of the product of operators A^ n A k is again exp ^ — ^-^J 



, the same as that in the corresponding term with 



%\k\ 

operators a^ n a k in Eq. (|4|). Intutively this is the correct way of representing the interaction of the ion and the laser 
in a q-analog harmonic oscillator trap. This form of the operator shows that the energy exchange of the center of 
mass motion of the ion occurs as the ion moves up or down in the trap. Since E n+ i — E n is n dependent, it implies 
an n dependent entanglement of the center of mass motion and the internal degrees of freedom of the two level atom. 
Using Eq. (|Io|), we can rewrite Eq. ( |l3|) as 

n=0fc=0 

The matrix elements of this operator are well defined(Eq. 14 of ref. jl[])and given for m < n by 

m * i*> = '^' fc) r '"'•' £; „, (£n y „ , m 

' ]J! £2*!(»-m + *)![n>-*],l 



,(*). 



Comparing F q with F , we notice that aside from the factor e V / , the effective lamb Dicke parameter for the 
loss and gain of motional state energy in a given interaction process in a g-analog trap is ef(N), that is it depends 
on the number of q-oscillator quanta linked with the state at a given moment. For the initial state of the system 
considered in our paper that is 

1 °° a™ 

\9,a) q = —j==Y^ —j= \9,n) , (16) 



expif 1 n=o \J[n] q \ 

the calculated value of q (a\ f 2 (N) \a) is 1.0004 , for a = 4 , e = 0.05 and r = 0.003. Of course it increases with 
increasing q. We take this opportunity to correct a typo in the 5 t,l line of Sec. V of our paper Q] and give again the 
parameters used for the numerical calculation. The values are e = 0.05, ZJ — £ = 50 and A = -4 = —50. In Eq. (|lT 



if we approximate ef(N) by e q — q (a\ f 2 (N) \a) q which for the case cited above gives e q — 0.05001, we can write 
the effective ion laser interaction operator as 



A 2 2 )e\ 2 / e ie "° ' e le « a = eV 2 



F g = eV ' * ) e \ * J e ^ a ' e ^ = e \ * - J e ^+ ») . (17) 
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We agree with the authors of the comment that the operator (A^ + A) is not a self adjoint operator and the basis 
vectors |m) are not analytic vectors of the operator exp (ie(A^ + A)^j. In our article |l| we start with a form for ion 
in ^-analog harmonic oscillator trap Hamiltonian that is analogous to the Hamiltonian of Eq. ([!]), but in the third 
paragraph of Sec. Ill of our article we make an approximation and put exp[ie(A^ + A)] = e^^ - ^' e ttA ^ e leA . We 
recall here that the Baker Hausdorff theorem (Eq.(^)) can not be applied here in case exact commutation relations 
for the operators A and A^as given by Eq.(^) are used. Only for the case q = 1, the operators A^ and A commute 
with L4t,A|. So the approximation amounts to using [AT, A] = — 1 at a certain level. The same approximation is 
used again to approximate e^^ - ^ by e^^~\ For the numerical calculation, the matrix elements of Eq. (|l5| ) 
go into the calculation of matrix elements of the type (e, m\ F q o~ + \g, n) . These points have been mentioned but not 
discussed in detail in our paper JjJ . The final form of the model Hamiltonian as used for numerical calculations is the 
one given by Eqs. ([n] and |lj) and not the one containing the operator exp([ze(v4^ + A)] . It is relevant to recall here 

that e\ 2 / e leA ^ e leA ^ e ^ 2 ' e leA e leAl due to Eq. (||). We conclude that the numerical calculations presented in 
ref. 0] are perfectly valid. 
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